Collective Excitations in High- Temperature Superconductors 
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Collective, low-energy excitations in quasi-two-dimensional d-wave superconductors are analyzed. 
While the long-range Coulomb interaction shifts the charge-density-wave and phase modes up to 
the plasma energy, the spin-density-wave excitation that arises due to a strong local electron-electron 
repulsion can propagate as a damped collective mode within the superconducting energy gap. It is 
suggested that these excitations are relevant to high-T c superconductors, close to the antiferromag- 
netic phase boundary, and may explain some of the exotic features of the experimentally observed 
spectral-density and neutron-scattering data. 
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An important aspect of superconductivity is concerned 
with collective modes as they may modify low-energy 
properties of superconductors Bj. In principle, these 
excitations are either attributed to a symmetry trans- 
formation under which the system is invariant or to a 
spontaneously broken continuous symmetry. One might 
then expect that, in addition to the usual charge and 
spin fluctuations, collective modes associated with the 
phase and the amplitude of the superconducting order 
parameter are important. Yet, in conventional supercon- 
ductors with a momentum-independent energy gap, col- 
lective modes have essentially no practical significance 
regarding the low-energy properties of the superconduc- 
tors mainly because the long-range Coulomb interaction 
causes the phase mode to appear at the plasma energy [|| . 
The only mode that is not diffusive at long wavelengths 
is the amplitude mode, but it has no direct coupling to 
charge or spin degrees of freedom making its observation 
demanding |jj. 

Moreover, excitonic states that are bound pairs of 
quasiparticles have been predicted to exist in supercon- 
ductors HI. Such states appear in angular- momentum 
channels other than the one in which the Cooper pair- 
ing occurs. Theoretically, they should be present once 
the effective electron-electron interaction has an attrac- 
tive partial-wave component with a given angular mo- 
mentum. Nonetheless, there is no experimental evidence 
for these kind of states possibly because of their small 
binding energy ||. 

While for most phenomena collective excitations in 
conventional superconductors can be ignored — for in- 
stance, the superconducting energy gap in the elec- 
tronic spectrum at the Fermi energy is not affected by 
these modes -- there are both theoretical and exper- 
imental reasons to expect that high-T c superconduc- 
tors may behave differently in this respect. Theoret- 
ically, the energy gap has a strong momentum depen- 
dence and there is a large local electron-electron repul- 
sion which may qualitatively change the nature of col- 
lective excitations and allow new ones to develop that 
are not related to the broken gauge symmetry. Ex- 



perimentally, neutron-scattering studies H in supercon- 
ducting Lai.6- a; Ndo.4Sr a; Cu04 have established the ex- 
istence of elastic peaks in the magnetic structure factor 
at wave vectors (| ± e, |) and (|, \ ± e) (measured in 
units of 2ir/a), providing direct evidence for (diffusive) 
spin-density-wave fluctuations in high-T c superconduct- 
ing materials. These excitations can be regarded as a 
manifestation of fluctuating charge stripes and antiphase 
spin domains which have shown to provide a natural 
explanation for the unusual features observed in angle- 
resolved photoemission experiments [Q . We may also ar- 
gue that they account for the broad "bosonic" feature in 
the electronic spectral density near the (^,0) and (0, i) 
points seen only in the underdoped materials Q. 

In this Note, we examine under what conditions col- 
lective excitations develop in d-wave superconductors and 
what are their experimental implications, when the quasi- 
particle picture is appropriate. Our most important ob- 
servation is that (i-wave superconductors close to the 
antiferromagnetic phase boundary support a low-energy 
spin-density-wave mode. As a consequence of the super- 
conducting energy gap, this mode can propagate coher- 
ently with reduced damping, unlike its precursor in the 
normal state where only a diffusive spin-density wave is 
realized. It can be excited by magnetic processes making 
it observable, for instance, by inelastic neutron scatter- 
ing. The mode is "massive" because rotational symme- 
try in spin space is not broken. Under specific conditions, 
however, the mass of the mode may vanish and may even 
become negative signaling an instability of the supercon- 
ducting ground state against a spontaneous creation of a 
spin-density-wave state. 

Consider, for example, the Hamiltonian 
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where n r = ^ CT ipl a ipra is the electron number operator 
at site r, t rr i is the tunneling-matrix element between 
sites r and r', fi is the chemical potential, and v(r) is 



an instantaneous electron-electron interaction. In de- 
scribing superconducting order, it is useful to express 
the Hamiltonian in the form H = Hbcs + -Hint, where 
the BCS and interaction Hamiltonians are 

#bcs = J2 *I( £ k f 3 - A k fi)* k (2a) 

k 
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Here, ^ k = (tp^ f/'l k i ) T is the Gor'kov-Nambu spinor, 
e k = J2r ^roe _qr — n is the single-particle energy rela- 
tive to the chemical potential, and v(q) = J2 r v(r)e~ q ' r . 
The fermion operators in real and momentum spaces 
are related by the unitary transformation tp ra = 
N~ 1 / 2 ^2 k tpk cr e lk ' r , where N is the number of sites in 
the system. Also, suppose that there exists a wave vec- 
tor Q such that £k+Q — — £k; «-e., the Fermi surface 
is approximately nested. This kind of situation may 
qualitatively arise in lightly doped high-T c superconduc- 
tors. For illustrative purposes, let only the nearest- 
neighbor tunneling matrix element be non-zero so that 
e k = — ^W (cos k x a + cosfc y a) — /i, where W is the half 
bandwidth and a is the lattice spacing. Therefore, at half 
filling (fi — 0), Q = Qo, where Q = (ir/a, ir/a). 

The energy gap A k is determined by requiring that 
the interaction Hamiltonian does not give any self-energy 
corrections to the energy gap || . This condition leads to 
the gap equation, 



A k = -^ U (k-p)Trf 1 G(p, 



(3) 



where G(p,u) = l/(wf — e p f 3 + A p fi) denotes the 
Green's function of the BCS Hamiltonian. We use the no- 
tation in which ^ = TV -1 5Z_ /(dw/27r) andp = (p, u)). 
In order to determine whether the system can support 
collective excitations, consider an effective two-particle 
interaction V that describes mutual scattering of two 
quasiparticles. The poles of the effective interaction then 
yield the energy and lifetime of two-particle collective 
excitations. In the ladder approximation, the Bcthc- 
Salpeter equation for T may be written formally as 



which is equivalent to the equation 

t(k, k'; q) = f (fc, k'; q) + i^2 T o(k,p; g)A(p, q)t(p, k'; q), 

p 

(4) 



with k(p, q) = G(p + q/2) ® G(p - q/2) T ; here, q is the 
total four-momentum of a quasiparticle pair. The bare 
interaction vertex, denoted as 



+ 



is T (k,k';q) = v(k - k')f 3 ® f 3 - v(q)P, with P = 
(fo — n) ® 0. Moreover, the outer product, [A (g) B] a i, = 
AijBki, where a = (ij) and b — (kl), is defined so 
that there is a one-to-one correspondence between the 
ordered index lists a,b £ {1,2,3,4} and (ij),(kl) 6 
{(11), (22), (12), (21)}. 

On a square lattice, it is convenient to define orthogo- 
nal functions ?7 Q (k) in terms of which the vertex func- 
tions and the electron-electron interaction can be ex- 
panded pi. For example, i>(k — k') = U + ^[cos^ — 
k' x )a + cos{k y — k')a] may be written as w(k — k') = 
J2 a VaVaityrjatk 1 ), where, for a = 0, v a = U is the on- 
site electron-electron interaction and, for a = 1,...,4, 
v a = V\/2 is the nearest- neighbor electron-electron inter- 
action. Below, we will assume that the on-site interaction 
is repulsive, U > 0, and the nearest-neighbor interaction 
is attractive, V\ < 0, as should be appropriate for a phe- 
nomenological model describing high-T c superconduc- 
tors. Similarly, t(k,k';q) = J2 a p 1^(9)77^)77,0 (k') and 

r (M';p) = Ea/jrS(«)^„(k)^(k'), with f<g(g) = 

Va^a/3T3 T3 — v(q)5 a o5poP. Thus, the Bethe-Salpeter 
equation becomes 

T aP (q) = T^(q) +X)fW(«)A rr ,(«)f y/9 («) ) (5) 



where A a/3 (q) = iJ2 P r la(p)Hp, q)rip (p)- To compactify 
the notation, define new matrices such that [r] Q/ 3 = T a p, 
etc. Then, it is immediately clear that the effective in- 
teraction vertex is 



f(< z ) = [l-f(°)( (Z )A( 9 )]- 1 f( )( (Z ), 



(6) 



where q — (q, f2 + i0 + ). Because, for q = 0, A(— p, q) — 
A(p,q), there is no mixing between even and odd par- 
ity sectors of T a p (q) : the subspaces {?7 k ,%. '} and 
{C k } describing scattering in the singlet and triplet 
channels decouple. In contrast, for q = Q, the subspaces 

i^k 'Ck } an d {Vk } are decoupled. When the sys- 
tem has particle-hole symmetry at the Fermi energy and 
fl ~ 0, a further factorization can be shown to occur; 
namely, particle-hole and particle-particle channels de- 
couple. This means that, for q = Q, the subspace {Q } 
decouples from {?7 k , Q '} and, for q = 0, the same oc- 
curs for {?7 k ~ } and {?7 k , ?7 k '}. This is particularly con- 
venient because 7? k determines the ci-wave gap function, 
A k = (A /2)r] k ' . At zero temperature, Eq. (||) becomes 



1 = -(VyA0£ k [^ -) ]74£k. Here, E k 



k . nere, £/ k = v4+^k is 
the quasiparticle energy in the superconductor. We find 
A ~ AWer 1 ^^^ , where W is the half bandwidth and 
Nf = N(fi) is the density of states at the Fermi energy 
in the normal state. 

The collective excitations are described by the poles of 
the effective two-particle interaction, Eq. (|6|), and they 
are conveniently classified by the symmetry properties 
of the system. The Hamiltonian (2) possesses a number 
of symmetries which lead to conserved currents. First, 
gauge symmetry, generated by the transformation ^ k — ► 
^k = elipT3 ^)i, yields charge conservation. On general 
grounds, one then expects that there is a collective mode 
associated with the phase of the superconducting order 
parameter. However, due to the long-range nature of 
the Coulomb interaction, the Anderson-Higgs mechanism 
shifts it at small momenta to the plasma energy. Second, 
the symmetry transformation, ^ k — ► ty' k = e^ 1 " 1 ^, is 
associated with an amplitude mode of the order param- 
eter pip . In contrast to s-wave superconductors, this 
mode is always over-damped in d-wave superconductors. 
Third, spin-rotational symmetry leads to a new mode in 
the particle-hole sector, which is driven by the on-site 
Coulomb repulsion. As an example, consider a cylindri- 
cal Fermi surface and a wave vector q = Q nesting two 
given k points with vanishing quasiparticle energies. For 
U ~ |Vi|, the energy Oq of the collective excitation is 
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when > 0. For A 0j |/x| « W, V' - w E k lk ~ 
IVi!" 1 . The parameter (3 = - [A 31os f (e) ] £ ^ mea- 
sures the magnitude of the particle- hole symmetry break- 
ing in the density of states N(e) at the Fermi energy 
and kq describes the contribution due to the mixing 
between the particle-hole and particle-particle channels. 
For U, \Vi\ -C W, it can be expanded as kq = «q +«q , 

where «g> ~ 2(%f \^f and «g> ~ (%)* (&)' ■ 
For /3 < 0, the collective excitation remains massive 
(f^Q > 0) for all values of Vq/U, down to the point where 
the superconductor becomes unstable. The excitation is 
damped, because it overlaps with the quasiparticle con- 
tinuum. For the other values of q, the continuum does 
not necessarily start from zero but at some finite value 
■dq = min k (_B k -|-i? k+q ). It is then possible to have collec- 
tive excitations with an infinite lifetime, when fi q < ?? q . 
The exact nature of the collective excitation is deduced 
by examining its coupling to the spin operator 5 z (q) — 

S k * k 7^ *k+q, where % — \t . That the mode is a 
spin-density wave becomes evident by computing the spin 
correlation function, x( c Lt t ) — {T T S z (q, t)S z (— q, 0)), 
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FIG. 1. The energy f2 q of the spin-density-wave excita- 
tion as a function of wave vector q at zero temperature for 
U/W = 0.485 and Vx/W = -0.608. Here, the superconduct- 
ing energy gap is Ao/W = 0.1 and the chemical potential is 
chosen so that the density of holes equals 10% (relative to 
half filling). The excitation energy Jl q is computed numeri- 
cally from the exact result, Eq. (ph. 



Specifically, its temporal Fourier transform is given by 
the formula 



x(«) = (7z|A(«)[l + r(«)A( g )]| 7 , 



(8) 



where the column vector \"/ z ) is defined as (aa\j z ) = 
^[fo]y<5 Q o- In the limit Q — > 0, our result reduces to 
the form obtained by non-conserving approximation p2[ , 
when the energy spectrum has particle-hole symmetry at 
the Fermi energy. However, for energies £1 > Ao or in 
the absence of particle-hole symmetry, one must use the 
general result, Eq. (g). Note that, by incorporating the 
mixing of the particle-hole and particle-particle degrees 
of freedom, kq accounts for, for example, the effect of 
any two-particle, spin-triplet excitations ]13[ ] at the mo- 
mentum Q. Similar calculation shows that no resonance 
develops for the charge response near q = Q. 

The most favorable conditions for observing these ex- 
citations are most likely found in undcrdoped high-T c 
superconductors close to the antiferromagnetic phase 
boundary. In the antiferromagnetic phase, the above col- 
lective mode is replaced by the Goldstone mode of the 
antifcrromagnet. Interestingly, for U > Vo, the system is 
unstable against a spontaneous creation of quasiparticle- 
quasihole virtual bound states at q = Q. This implies 
a phase transition to an antiferromagnetic state. In con- 
trast to d-wave superconductors, the Ward identity [p|Jll|l 
excludes the spin-density- wave collective mode in conven- 
tional superconductors with a momentum independent 
gap function. 

To obtain a quantitative understanding of the disper- 
sion relation of the collective mode, we resort to numer- 
ical methods. Figure 1 shows the energy of the spin- 
density wave as a function of wave vector near half filling 
with hole density equal to 10%. The minimum value 
of the excitation energy, £l q /2Ao ~ 0.08, is located at 
the wave vector q ~ (0.957r/a,7r/a), implying that the 



superconductor becomes unstable against a spontaneous 
creation of an antiferromagnetically ordered state with 
vertical (horizontal) antiphase domain walls. In con- 
trast to neutron scattering data in superconducting 
Lai.6- a; Ndo.4Sr a ;Cu04 which show that the mean separa- 
tion £ between antiphase domain walls should scale with 
the hole density ias!~ a/2x, we find that £ is more than 
by a factor of two longer than the experimental one. (In 
the present approximation, £ also depends on Ao.) The 
failure to predict correctly the domain-wall periodicity is 
similar to the problem of describing the static, incom- 
mensurate stripe order in the Hartree-Fock approxima- 
tion ]ljj . Although at zero temperature exciting vertical 
spin-density-wave fluctuations requires the least amount 
of energy, the rotation of their orientation relative to the 
underlying lattice constitutes a relatively soft mode; see 
Fig. 1. For example, the excitation energy at the saddle 
point q = 0.98Qo is about 60% larger than the minimum 
energy required to excite a vertical mode. The form of the 
dispersion relation £l q is affected by the lifetime effects: 
the troughs clearly visible in Fig. 1 mark the boundaries 
between damped (f2 q > t? q ) and undamped (f2 q < t? q ) 
excitations. For example, the excitations with q = Qo 
and 0.94Qo are undamped whereas the excitation with 
q = 0.97Qo is damped. 

Figure 2 illustrates the behavior of the imaginary part 
of the spin correlation function, Eq. (ra). The collective 
excitation produces a distinctive resonance structure in 
x"(q, fi) at energies specified by the dispersion relation, 
fi = Sl q . At low energies, the resonance becomes nar- 
rower as the excitation energy f2 q decreases. However, 
in the vicinity of the (it /a, it /a) point, for example, the 
collective excitation appears below the onset energy of 
the two-particle continuum $ q . Accordingly, the collec- 
tive excitation would acquire an infinite lifetime yield- 
ing a resolution- limited peak in x"(q, $7) — also, below 
this onset energy, the usual quasiparticle contribution to 
x"(q, fi) would vanish — if the lifetime of quasiparticles 
were infinite. Clearly the collective mode is an important 
new feature describing the low-energy spin correlation 
function |15| . 

In a d-wave superconductor, strongly-scattering impu- 
rities induce virtual bound states Pq| by modifying the 
local potential energy of electrons at the impurity site. 
In addition, they may also change local Coulomb in- 
teraction between electrons at the impurity site. Such 
an effect can be accounted for by including the term 
t^imp'raroT ri ro!, where r = ro is the location of the im- 
purity. A straightforward calculation shows that this 
effect leads to a virtual bound state with the energy 
ft /2A ~ 2 v /£/ c /£/ imp - 1, where U c ~ 1/ttN f . In con- 
trast to an impurity potential, the impurity interaction 
has a critical value U c , below which the ground state is 
nonmagnetic and the spin quantum number of the reso- 
nance state equals to ^i an d above which the impurity 
becomes magnetic in the sense that one electron spin is 
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FIG. 2. The imaginary part of the spin correlation func- 
tion x"(qi^) as a function of (a) the wave vector q with 
Q/2Ao = 0.13 and (b) Q and q x with q y = n/a at zero tem- 
perature for U/W = 0.485 and Vi/W = -0.608. A small 
broadening of resonances is obtained by surmising a finite en- 
ergy resolution of magnitude ■y/W = 3 ■ 10~ 3 due to quasipar- 
ticle lifetime effects. Here, the superconducting energy gap is 
Ao/W = 0.1 and the chemical potential is chosen so that the 
density of holes equals 10% (relative to half filling) . 



trapped to the impurity site. 

Finally, one may ask whether excitonic states of bound 
pairs of quasiparticles are feasible in ci-wave superconduc- 
tors. It is immediately clear that while the effective in- 
teraction has an attractive partial wave in the (extended) 
s-wave channel, it is nevertheless too weak near half filling 
to support any excitons with the same angular quantum 
number. It is only far away from half filling that these 
states might appear as virtual bound states because of 
the proximity to a superconducting state with extended 
s-wave symmetry. 

In conclusion, we have shown that d-wave supercon- 
ductors can support propagating collective modes that 
are best described as spin-density waves. It is then nat- 
ural to anticipate that the most favorable conditions for 
detecting them are found close to the antifcrromagnetic 
phase boundary. Furthermore, superconducting fluctu- 
ations couple particle-hole and particle-particle excita- 
tions allowing the latter ones to be probed by the usual 
means. This coupling is particularly important if the two 
excitations are (nearly) degenerate. 
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